
######绘制生存曲线
setwd("E:\\1.Ovarian_Treg\\1.picture\\6.survival")
list.files()
#install.packages("survminer")
s2.exp=read.table("rs_tcga.txt",sep = "\t",header = T,row.names = 1,check.names = F)
library(survminer)
library(survival)
os.time2=s2.exp$futime
status2=s2.exp$fustat
############################################
#?ggsurvplot
res.cut=surv_cutpoint(s2.exp,time = "futime",
                      event = "fustat",variables = "riskScore")
res.cat=surv_categorize(res.cut)
fit=survfit(Surv(os.time2, status2)~riskScore, data=res.cat)
col=ggsci::pal_npg("nrc")(10)
ggsurvplot(fit, 
           surv.median.line = "hv", # 添加中位数生存时间线
           
           # Change legends: title & labels
           legend.title = "TCGA-Risk", # 设置图例标题
           legend.labs = c("high", "low"), # 指定图例分组标签
           
           
           # Add p-value and tervals
           pval = TRUE, # 设置添加P值
           pval.method = TRUE, #设置添加P值计算方法
           conf.int = F, # 设置添加置信区间
           
           # Add risk table
           #risk.table = TRUE, # 设置添加风险因子表
           #tables.height = 0.2, # 设置风险表的高度
           #tables.theme = theme_dark(), # 设置风险表的主题
           
           # Color palettes. Use custom color: c("#E7B800", "#2E9FDF"),
           # or brewer color (e.g.: "Dark2"), or ggsci color (e.g.: "jco")
           palette = col[1:2], # 设置颜色画板
           # ggtheme = theme(panel.border = element_rect(fill=NA,color="black", size=1, linetype="solid"))+
           # 
           #   theme(panel.grid = element_blank()),# Change ggplot2 theme
           title = "Survival curve", # 添加标题
           font.main = c(16, "bold.italic", "black"), # 设置标题字体大小、格式和颜色
           font.x = c(14, "bold.italic", "black"), # 设置x轴字体大小、格式和颜色
           font.y = c(14, "bold.italic", "black"), # 设置y轴字体大小、格式和颜色
           font.tickslab = c(12, "plain", "black"),# 设置坐标轴刻度字体大小、格式和颜色
           linetype = "solid",#改变线条的渐变
           legend=c(0.8,1))##图例位置


ggsave("TCGA.pdf",height = 3.2,width = 3.5)
